Model of interacting dark energy in five-dimensional 

brane cosmology 



Mubasher Jamil 

Center for Advanced Mathematics and Physics 
National University of Sciences and Technology 
E&ME campus, Peshawar road, Rawalpindi - 46000, Pakistan 

January 16, 2009 



Abstract 

We here investigate the model of interacting dark energy in the context of five dimen- 
sional brane cosmology. The effective equations of state of dark energy are evaluated for 
various choices of the variable time dependent cosmological constant. We have found that 
the interacting dark energy obeys the phantom divide/crossing scenario in this generalized 
model. It is also shown that interacting dark energy in this generalized model also resolves 
the cosmic coincidence problem. 

Keywords: Dark energy; Chaplygin gas; cosmic coincidence problem; cosmological constant; 
dark matter; phantom energy; phantom crossing. 



1 Introduction 



Recent astrophysical observations give bunch of convincing evidence of our universe under- 
going accelerated expansion [HEJEIIHEIEIEIIE]- ^ perturbing feature of this phenomenon 
is that it was preceded by a decelerated expansion, so we need to ask what caused this 
sudden transition [9]. Observations also show that this shift is rather marginally recent 
(less then one Gyr), hence it poses why it happened so recently. If we assume that the 
source which is driving this expansion is some mysterious 'dark energy' then one needs to 
ask what is the composition of this exotic matter, also why it has become dominating all 
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of a sudden at present time (see [10|, [XT], [12], [T3l [H] for reviews on dark energy). In other 
words, why it was negligibly small in earlier deceleration phase. Investigations of WMAP 
show that energy densities of dark energy and matter are almost comparable at present 
time. This leads to a problem named as the 'cosmic coincidence problem' (CCP). The 
problem aggravated when it was shown that radiation energy density was also equivalent 
to that of dark energy, leading to a 'cosmic triple coincidence problem' |15j . however it is 
recently addressed in the context of triple interacting fluids [15]. Further if the universe is 
dark energy dominated, will it expand forever or may decelerate at some instant as well. 

The dark energy is generally represented by a phenomenologically motivated equation of 
state (EoS) p de = u de p de , where p de and p de are the pressure and the energy density of the 
dark energy, interlinked by a dimensionless parameter u de [17]. In order to produce the 
accelerated expansion, it requires u de < —1/3. In recent years, several theoretical models 
have been proposed to understand the nature and dynamics of dark energy, however 
almost all these models either require fine tuning of their model parameters or yield 
quantum or gravitational instabilities that are needed to be removed. Most prominent 
dark energy proposals are based on cosmological constant A [TBI EE], quintessence (2TJJ, 
[21] . k-essence [22], [23], phantom energy |24j, quintom model (25] [26], geometric dark 
energy [27], holographic dark energy [28] and tachyons [29] ED], to name a few. It has 
been pointed out that quantum effects can yield a super-accelerated phase of cosmic 
expansion (without any need of introducing ghosts, phantoms or tachyons) and that these 
quantum effects yield stable solutions [31], [32], |33] . The precise determination of u de is a 
more challenging and interesting problem in itself. Recent observational data gives the 
estimate —1.67 < uj de < —1.05 at 95% confidence level [33]. It also supports the notion 
of an evolving u de , hence it requires the parametric form u de (z), where z is the redshift 
parameter [35], [36] [37J [38] [391 SO] ■ For larger redsifts (quintessence dominated), u de > — 1 
while at some instant u de = —1 (A dominated era) and later it is u de < — 1 (phantom 
regime). Therefore the CCP is rephrased as 'why now u de = — 1?' 

In the last few years, the CCP is addressed by invoking a non-minimal interaction 
between dark energy and dark matter (or simply 'matter' for convenience), the so-called 
interacting dark energy model @Il[l2l[I3[32[l5l[^ The model 

involves a coupling parameter and an energy exchange term to govern the interaction. 
The energy exchange term is adjusted so as to satisfy the global conservation law for 
the interacting system. Moreover the interaction is dynamic i.e. the energy is exchanged 
between the interacting components with equal degree of freedom. It is exactly this feature 
that helps in maintaining the equilibrium of densities of the interacting components i.e. 
the ratio of energy densities roams around to unity. It is recently shown that the coupling 
parameter cannot be negative in order to avoid possible violation of the cherished second 
law of thermodynamics, however, small positive values are permissible to account the 
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decay of dark energy into matter [5l[ |55[ 56J. The notion of decay of dark energy into 
matter particles is well motivated from the theoretical arguments and henceforth predicts a 
matter dominated universe. Hence a universe governed by the interacting components can 
undergo deceleration phase preceded by an acceleration phase. Further, the interaction 
also saves the universe from undergoing an imminent 'big rip' (tearing apart of spacetime 
structure leading to a future spacelike singularity). This model also favors a bouncing 
universe (free from cosmological singularities like big bang, big rip, big crunch etc) since 
the model forces the components to interact, thereby avoiding cosmological over-densities. 

We here extend our earlier work on interacting dark energy [53, [58] in the context of 
five dimensional brane- world model [59]. Last few decades have seen a considerable ad- 
vancement on the theories of extra dimensions and have revealed deep insights about the 
structure of spacetime, elementary particles and forces of nature [60]. It has been sug- 
gested that the mysterious dark energy is also a manifestation of extra spatial dimensions 
|61j . We investigate the behavior of dark energy in the brane world model and determine 
various EoSs for the dark energy for different choices of the time dependent cosmologi- 
cal constant. All the effective EoS describe the phantom crossing scenario under certain 
conditions. It is also discussed that this generalized interacting dark energy model fairly 
resolves the CCP. 

2 Interacting dark energy model 

We start by assuming the background to be spatially homogeneous and isotropic Friedmann- 
Robertson- Walker (FRW) spacetime, in the context of brane-world gravitation model, 
which is given by [62} [59] 



Here A; is a curvature parameter which refers to a spatially spherical (k = +1), Minkowskian 
(k = 0) or hyperbolic (k = —1) spacetime. Note that these cosmological models are also 
correspond to closed, flat or open respectively. Also a(t) is the dimensionless scale factor 
while ( is an arbitrary function of time t and we choose it £(t) = a n , whereas n is a 
constant. Moreover ip is the fourth spatial dimension. The spacetime is further assumed 
to contain two fluids namely matter and dark energy. The corresponding energy den- 
sities are p m and pd e while the respective pressures are p m = (pressureless dust) and 
Pde = ^dePde 7^ 0. The combined matter energy distribution is given by a perfect fluid 
stress energy tensor 



ds 2 = dt 2 - a 2 (t) 



dr 2 





1 — kr 2 




(2) 
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Here is the five-velocity vector which satisfies u^u 11 
first FRW equation is 



(n + l)H 2 + h- 



1 



3M p 2 



(Pde + Pm) + 



1 with = 0,1,2,3,4 . The 
A(t) 



(3) 



Here M p = (87rG) 1 is the reduced Planck mass. The corresponding energy conservation 
equation is 



Pde + Pm + (3 + n)\p m + Pde (l + oo de )}H = -M 2 A, (4) 

which can be decomposed into two non- conserving equations for both matter and dark 
energy as 



Pde 



(3 + n)(l + u de )p de H 
(3 + n)p m H 



-Mi 



Q, 



,A 



Pr, 



-m 2 ;- + q. 



(5) 
(6) 



Note that addition of the above two equations leads to the energy conservation (4). Above 
Q is the energy exchange term for the interaction. We here choose Q = 3Hb(p m + p de ) 
[63], [64]. where b is the coupling parameter (or transfer strength). Due to unknown nature 
of both dark energy and dark matter, the interaction term can not be derived from the 
first principles. It is worthy to note that if Q < than it will yield the energy density 
of dark energy to be negative at sufficiently early times, consequently the second law of 
thermodynamics can be violated [65] hence Q must be positive and small. Because of 
the underlying interaction, the beginning of the accelerated expansion is shifted to higher 
redshifts. 



Further the density parameters are 



n rae 
1 ide — 



Pr, 

Pc 
Pde 

Per 
k 



Pde 



0; 



'H 2 ' 



Using the above parameters in Eq. (3), we obtain 



(/(:' 



1 + n + flu — 



A(t) 
3H 2 ' 



We here define the dimensionless ratio of densities 

Pm f^rr 



Pde 



n 



tie 



(7) 
(8) 

(9) 
(10) 
(11) 
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To check how this density ratio evolves with time, we differentiate it w.r.t t to get 



Pr, 



Pr. 



Pde 



dt Pde iPm 

Using Eqs. (5) and (6) in (12), we obtain 

T 1 + r 



= f(r x ). 



(12) 



r x (3 + n)H 



UJde + 



(3 + n)H r a 



+ 



A(r, - 1) 
QH 3 Vt m {3 + n) 



where 



T = 3Hb{l + r x ) 



(13) 



(14) 



is the decay rate, related to Q — Tp de . Thus as r x approaches 1, the last term on right 
hand side in equation (13) becomes negligible. It is also termed as the 'soft coincidence' 
since \r x /r x \ < H [66]. Further if n = then Eq. (13) reduces to the one discussed in 
J. Using Eq. (14) in (13), we obtain 



r x {3 + n)H 



36(1 + r x 



+ 



A(r a 



r x (3 + n) QH 3 n m (3 + n) 
The critical points (or stationary solutions) are obtained by solving r x = to get 



(15) 



A + 18bH 3 n r 



+ r x 



6H 3 Q m {6b + Ude(3 + n)} - A 



+ 18bH*n r 



0. 



(16) 



The above equation yields two roots as 

1 



r x ± 



[A - QH s n m {Qb + uj de (3 + n)} 



2 (A + 18bH 3 n r 

±J -72H 3 Q m b{A + 18bH 3 Q m ) + {A - 6H 3 [6b + (3 + n)w de ]fi m } 2 ]. (17) 



It is recently shown in [50J that any model of interacting dark energy can resolve the 
cosmic coincidence problem if the function f(r x ) satisfies 



df 

-(r = r xi )<0, 



;is) 



where r xi for i = 1,2, ... are the roots of f{r x ) = 0. It needs to be stressed that not all 
roots will satisfy (18) but those which do satisfy it, are termed 'stable equilibrium points'. 
In our model, the condition (18) together with (17) yields 



f(r x± ) 



-72H 3 Q m b(A + 18bH 3 Q m ) + [A - 6H 3 {6b + (3 + n)u de }fl r , 



6H 2 tt r 



(19) 
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Hence it is clear that f'(r = r x _) < and r x _ is the only stable equilibrium point of 
our model. Thus interacting dark energy model in brane-world gravitation theory fairly 
alleviates the CCP. Our next task is now to determine the effective EoS for dark energy. 



The parameter r x in (11) is related to the density parameters (7) - (9) as 



1 

&>de 



1 + n + Q k - Q de - 



A(t) 

m 2 



We further define the effective equations of state for dark energy and matter as 

1 r 



eff 
de 



r 

3H~ 



^de + 7TT7) U m 



eff 



From Eq. (5) we have 



^de 



-1 - 



(3 + n)Hp de 



Q 



r r 3H' 



M 2 — 
p 2 



(20) 



(21) 



(22) 



Using Eq. (22) in (21), we get 



00 



eff 
de 



-1 - 



or we can write 



eff 
de 



(3 + n)Hp de 



Pde 



A 



+ 



nT 



3(n + 3)#' 



(23) 



+ 



nb(l + r x ) 



(3 + n)Hp de 6(3 + n)Hm de ' n + 3 
We represent the dark energy by the modified Chaplygin gas (MCG) equation of state 



(24) 



Pde 



B 

Ap de - — 



(25) 



where A, B and a are constant parameters. The MCG best fits with the 3— year WMAP 
and the SDSS data with the choice of parameters A = —0.085 and a = 1.724 which 
are improved constraints than the previous ones —0.35 < A < 0.025 [68]. Recently it 
is shown that the dynamical attractor for the MCG exists at u de = —1, hence MCG 
crosses this value from either side u de > — 1 or u de < —1, independent to the choice 
of model parameters [69]. A generalization of MCG is suggested in [70] by considering 
B = B(a) = B a a , where o and B are constants. The MCG is the generalization 
of generalized Chaplygin gas p de = —B/p^ e [7U [72] with the addition of a barotropic 
term. This special form also appears to be consistent with the WMAP 5— year data 
and henceforth the support the unified model with dark energy and matter based on 
generalized Chaplygin gas [731 [74]. in the cosmological context, the Chaplygin gas was first 
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suggested as an alternative to quintessence and demonstrated an increasing A behavior 
for the evolution of the universe [75]. Recent supernovae data also favors the two-fluid 
cosmological model with Chaplygin gas and matter [76]. The density evolution of MCG 
is given by 

p de = (X + C 1 a Y )^, (26) 

where X = B/(l + A), Y = 3(1 + a)(l + A) and C\ is the constant of integration. The 
time derivative of p de is given by 

p de = -3(1 + A)C 1 #(X + C 1 a y )^a y . (27) 

Using Eqs. (24) to (27), we obtain the effective EoS for dark energy as 

e ff 3(1 + A) f _ X \ _ A nb(l + r x ) 

de ~ + 3 + n [ L f%") 6(3 + n)Hm de + n + 3 ' W 

Models with variable A(t) are physically more appealing and theoretically rich in pre- 
dictions as compared to constant A. If A < 0, it gives a decreasing behavior of A with 
time. Physically it may explain inflationary expansion at earlier times while an acceler- 
ated expansion in current time. Similarly, A > represents an increasing A with time. It 
can be best interpreted in a model of bouncing cosmology where a universe is free from 
any potential cosmological singularities like the big bang one and bounces back near the 
imminent singularity. In such a scenario, a smaller A corresponds to a deceleration phase 
followed by a smoothly evolving larger A which results in a de Sitter like expansion. If we 
assume that this later expansion is driven by an exotic phantom energy (u de < — 1), then 
the later one decays into matter particles creating a matter dominated universe again [77] . 
Thus in a bouncing universe, an otherwise big rip singularity is replaced by the matter 
creation scenario. Hence if u de < — 1, it eventually leads to two interesting results: first, 
the existence of a bouncing universe and second, the decay of dark energy into matter or 
the model of interacting dark energy [7H1 05] . It is recently suggested using inhomogeneous 
EoS for dark energy that the dark energy dilution becomes faster in de Sitter expansion 
which involves strong interaction between dark energy and matter [79]. We shall now pro- 
ceed to determine ufJJ by assuming dependencies of A on various cosmological parameters 
and determine conditions under which it will become super-negative. 

From Eq. (28) we see that u e JJ < — 1 if A > and n > —3. The scenario of dark 
energy dilution into matter arises for b > or b — > 1 if we restrict < b < 1 while 
a b = corresponds to a non-interacting dark energy model. It implies that for some 
specific values like n = — 1, —2 the last term in Eq. (28) will also be negative. The case 
for b < refers to matter decay into dark energy, which is not relevant here. Let us choose 
A(t) = C 2 tP than Eq. (28) yields 

eff _ 3(1 + 4)/^ X\ C 2 (3t 13 - 1 , nb(l + r x ) 



^ ~ 1+ 3 + n ^ pT) 6(3 + n)M de + n + 3 " ( 29 ) 
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Next we choose A(t) = C^e 7 * which gives 



eff_ , , 3(1 + A) ^ X ^ C 37 e 7i , nb(l + r x ) 



~ 1+ 3 + n ^ p^J 6(3 + n)M de + n + 3 " (30) 

Thus in the above two Eqs. (29) and (30), A > translates into (3 > 1 and 7 > 0. Next 
we take A(t) = C^a 5 which enables us to write 



«4" = -1 + 



Pde 

which can alternatively be written as 



3(1 + A) f X \ C 4 5a s nb(l+r x ) 

3 + n V Pat") 6(3 + n)M de + n + 3 ' 1 ' 



r 1 & i Y 

3(1 + A) ( X\ C ^ \h (Pdt a ~ X )\ nb{l+r x ) 



+ v x ' . (32) 



3 + n V P*/ 6(3 + n)# 2 n de n + 3 

Here we require 5 > and <5 > Y to get a super- negative EoS. If we take A(t) = C 5 H V 
then we have 

^f = . 1 + ^( l _4\.£^ + ^±lA. (33) 



3 + n V Pdt a J 6(3 + n)Q de n + 3 

The phantom crossing scenario is more prominently observed from H > (u de < — 1), 
H = (uo de = — 1) and H < (cj^ > —1). More specifically, at the transition uj de = 
u e JJ = —1, we require 



3(1 + A) f X 
n(l + r x )\ p d f 



6 = -3^(l-^)- (34) 



Since < b < 1 and n > — 3 for A > 0, we obtain a restriction — 3 < n < from Eq. (34). 
Lastly we take A(t) = C%p de hence we get 



ff = 3(1 + A) f _ _ C&p^P* nb(l+r x ) 

de 3 + n V pltV 6(3 + n)M de n + 3 



which we can be simplified to yield 

e//_ 1 + 6ggPL \ / X \ n6(l+r a ) 

-- 1 + iT^l 3 + 2F^JI v 1 -^J + ^T^- (36) 

At the epoch of phantom crossing, we require u de = — 1 = oo e JJ to get 

" = -snrb ( 3 + wk) " w) (37> 

Note that Eq. (37) is reduced to (34) if e = 0. 
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3 Conclusion and discussion 

In the present work we attempted to resolve the cosmic coincidence problem in the con- 
text of five- dimensional brane world gravitation theory. The CCP is fairly alleviated since 
stable stationary solution exists for the dynamical system. The dark energy is represented 
by the modified Chaplygin gas and is further assumed to interact with the matter. This 
interaction leads to a phantom crossing scenario. We have also determined various ef- 
fective EoS for dark energy using different choices of A(t), since EoS of dark energy will 
change if it interacts with matter. This paper also presents a generalization of the work 
in [80] where it is shown that the dark energy with the MCG EoS crosses the phantom 
divide in the background of four dimensional FRW spacetime. 

Models of interacting dark energy have taken considerable interest in recent years. The 
model mimics ACDM at early times while it gives a finite dark energy-dark matter ratio 
at late times so that coincidence problem is alleviated. It turns out that dark energy 
domination is merely a transient event and will be replaced by the dark matter dominant 
era once again [81J. It is recently shown that this model can resolve the cosmic age 
problem as well since simple dark energy cannot remove the problem [82] ■ The cosmic 
age predicted by the interacting model is predicted to be greater than the ACDM model 
which is consistent with the observations and hence alleviates the cosmic age problem. 
Further, this model is also investigated in the context of loop quantum cosmology a 
theory in which all cosmological singularities are avoided due to quantum effects [83]. It 
is also of some interest that the coupling parameter for the interaction yields a variable 
Newton's gravitational constant [M]. From the observational point of view, a small but 
non-vanishing interaction is reported from the analysis of the dynamics of 33 relaxed 
galaxy clusters like Abell A586 [85] . 
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